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1. Introduction 

The precise calculation of the charm and bottom spectra is an important goal of lattice QCD 
for several reasons: 

• There are many 'gold-plated' states: narrow, stable, and experimentally accessible. 

• The splittings in heavionium have particularly good properties for determining the lattice 
scale. 

• It is an important test of the actions used for b and c, which can then be used to calculate 
decay constants which in turn are crucial for determining CKM matrix elements. 

• Control of systematic errors are well developed in calculations of the B and D spectrum on 
the lattice due to relative insensitivity to heavy-quark polarisation effects and to light quark 
masses. 

The Highly Improved Staggered Quark (HISQ) action [|l|, allows unprecedented control of 
discretization errors in numerical lattice calculations. We use HISQ s and c valence quarks with 
NRQCD b quarks on MILC lattices with Nf = 2 + \ flavors of ASQTAD sea quarks to calculate 
the masses of the Bs and Be mesons. 

2. Heavy-light 2-point functions 

For increased statistics, we use random sources r]{x,ta), defined as a three-component random 
complex unit- vector defined on each point in the source time slice t^. These are the sources for 
the inversion of the HISQ strange and charmed valence quark propagators. The HISQ propagators, 
being staggered, are spinless and we need to convert them to 4-component propagators to combine 
with NRQCD propagators in a Z7-light correlator. This is readily done at the sink by multiplying by 
the standard staggered-to-naive transformation Q.{x) = Hi Tf ■ Since the propagator source disap- 
pears once the propagator is made, the Q. factors needed at the source (when a random wall is used) 
must be transferred to the source of the heavy quark propagator and how to do this is described 
below. 

NRQCD propagators are then made from a source which includes the same random wall with 
which the HISQ propagators are made, plus the Q. factors and in addition different Gaussian smear- 
ing factors of varying radii r, chosen to allow improved overlap with the ground state B^ and Be 
mesons so that their energies can be extracted accurately at early correlator times. We therefore 
initialise A^smear NRQCD b propagators by setting 

Gi{x,t =Q) = Y^S{\x-x\\ri)^{x)Q.{x) (2.1) 
x! 

and evolve with 

„ ^ (l - f ) (l - ^)"./(.) (.-§)"('-?) ^'M. (2^2, 
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We use an improved lattice NRQCD Hamiltonian |Q] : 

g ~ a^A^ a(5(2))2 

Finally, at each time-slice, we combine the A^smear NRQCD propagators and the valence s or 
c quark propagator with the same smearing functions at the sink end, giving us A^smear x A^smear 
and Be meson correlators. 

3. 2-Point function effective masses and noise 

The expression for the variance of the Bs correlator 

[{Gb, (/, j;t- ?o)Gb, (/, i; t-to))- {Gb, (/, j; t-to)f] (3.1) 

contains in it propagators for bsbs four-quark states. The lightest combination on the lattice is 
Vb + Vs- Therefore the error on the propagator falls like ^^2(^1?/,+^'?!)' while the signal falls like 
^-Mg.t ^ as can be seen in Figure [I[ So the signal-to-noise ratio degrades exponentially by the 
mass difference seen in the figure. This degradation is worse than for states and necessitates 
smeared sources to be able to fit propagators to small t values. Figure |l] also illustrates that effective 
masses of propagators GB,{i,i',t — to) with Gaussian smearing with radius 2 and 4 {G2G2 and 
G4G4) approach the Mb^ plateau rapidly compared to the local source-sink combination (LOLO). 

4. Extracting Mb, and Mb, 

We fit the measured Bs and B^ correlators to the form 

k=i k'=i 

where / and j respectively index the source and sink smearing functions. The second term is an 
oscillating parity partner state. We perform simultaneous Bayesian fits of the Gmeson(?,j;f — to), 
looking for stability with respect to fit range and A'exp- 

Since the NRQCD Hamiltonian does not include a mass term, there is a shift in the energy of 
the p = states relative to the continuum mass. To correct for the energy shift in Mb^ and Mb^. we 
use the relationship: 

Mb,^ = {eb,^-\e,j)^^^+Im,t;, (I) 

where Eb^ or Eb^ is the ground-state Ei . Mg^ on the right-hand side is the spin-averaged experi- 
mental masses of bb states: Mg^ = (3Mx +M^^)/4, where we use the recent BaBar measurement 
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Figure 1: Plot of effective masses of (local-local) Bs correlator and Bs correlator error. While the effective 
mass of the correlator matches the experimental Mb, (corrected for the energy shift), the effective mass of 
the correlator error gives j {M-qi^ +M,jJ. Several source-sink smearing combinations are shown. 

of the T(15')-T]fo(15) hyperfine splitting [^. is the corresponding spin-averaged lattice energy, 
calculated with NRQCD b quark propagators on the same configurations [^. For Mb^ we also 
explore two other methods for cancelling the energy shift: 

Mb, = [eb, - \iEtB + Ec-,))^+\ {M,-, + M,,) , (II) 

where Mcr = (3M,^ + Mf,J /4, and Ecr is the corresponding spin- averaged lattice energy. The final 
method to extract Mb^ is 

Mb, = {Eb, - {Eb^ + Ed, - + {Mb, + Md, - M^J , (III) 

where M^, (Erj,) is the mass (lattice energy) of a fictional ^'s pseudoscalar: M,j^ = ^2M|- — M|. 

We show results for several ensembles at three lattice spacings and different light sea quark 
masses for Mb, in Figure ^ and Mb, in Figure ||. The statistical error (shown in figures) is dominated 
by the uncertainty ina^^. We use MILC ri/a values [^, to set the scale ensemble-by-ensemble. 
We then convert to physical units with ri = 0.321fm [^]. An additional systematic error due to 
the 1.5% uncertainty in r\ must be added in at the end with other systematics. Expressions (||) 
and (p^) reduce the uncertainty from errors in (ri/a) by using {Eb, — ^{Ef^-j^ + Ecc))Mt and {Eb, — 
{Eb, +Ed, —Eri,))\att, which are small relative to the {)^^^^ quantity in (|). 

5. Conclusions 

While preliminary, these results show the potential for precise calculations of Mb, and Mb, 
using HISQ valence quarks and NRQCD b quarks on a 2 + 1 flavor ASQTAD sea. The lattice 
results show little sensitivity to sea quark mass. Method (|) and to a lesser extent (||) for Mb, 
show a small dependence on the discretization scale. A slightly mis-tuned in the coarse and fine 
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Figure 2: Lattice calculations forMg, on very coarse 16^ x 48, a ' « 1 .3 GeV, coarse 20^ x 64 and 28^ x 64, 
fl"' « 1.6 GeV and fine 28^ x 96, a"' « 2.3 GeV configurations. 
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Figure 3: Lattice calculations and experiment for Mb^ on very coarse, coarse and fine configurations. For 
clarity we do not show the calculations from method (I), with its significantly larger error bars. 



ensembles ||8|] may contribute some systematic eiror in Mg^ values derived from (|I|), and will be 
corrected for in future work. Relativistic corrections and electromagnetic effects are two possible 
sources of systematic error. Future efforts will improve statistics and include finer lattices. 
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